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INTERACTION  BETWEEN  SMALL  PERTURBATIONS  WITH  DISCONTINUITIES 
AND  THE  STABILITY  OF  SHOCK  WAVES  IN  MAGNETOHYDRODYNAMICS 

V.  M.  Kontorovich,  Khar'kov 
1.  Annotation 

We  will  examine  the  Incidence  of  a  plane  monochromatic  wave  on  a 
surface  of  discontinuity  in  magnetohydrodynamics.  A  simple  geometric 
method  is  proposed  for  plotting  divergent  (reflected  and  refracted) 
waves  which  occur  in  this  case  (V.  M.  Kontorovich  [1]). 

Based  on  the  results  obtained,  we  will  examine  the  problem  of 
the  resistance  of  shock  waves  to  splitting  (A.  I.  Akhiyezer  [2])  with 
respect  to  oblique  perturbations.  It  Is  shown  that  the  number  of 
waves  diverging  from  the  discontinuity  does  not  depend  on  the  angle 
of  incidence.  We  stress  the  need  to  classify  waves  into  incident  and 
divergent  waves  by  group  velocity  and  not  by  phase  velocity  since  at 
certain  angles  of  incidence  their  projections  onto  the  normal  to  the 
discontinuity  have  different  signs.  With  a  classification  based  on 
phase  velocity  we  would  have  the  paradox  of  unstable  shock  waves  in — 
magnetohydrodynamics. 

The  geometric  method  of  plotting  waves  diverging  from  the  dis¬ 
continuity  can  be  applied  in  simplified  form  to  the  interface  otf!  two.  t 
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media  in  a  compressible  fluid,  and  also  to  a  tangential  discontinuity. 

We  examined  the  occurrence  of  surface  waves  in  a  compressible 
fluid  upon  interaction  between  perturbations  and  shock  waves.  ,  It 
is  shown  in  particular  that  surface  waves  arising  upon  reflection  of, 
waves  from  the  interface  and  on  their  passage  through  the  boundary  in 
a  compressible  fluid  ( R.  N.  Roberts  [2])  are  a  limiting  case  of  that 
branch  of  magneto-acoustic  waves  which  approach  ordinary  sound  and 
disappear  in  an  infinite  fluid  when  S-» 

We  observed  the  peculiarities  of  the  Doppler  effect  during 
interaction  between  small  perturbations  and  shock  waves.  We  found 
frequency  drifts  of  the  arising  perturbations  relative  to  the  frequency 
of  the  Incident  perturbation  ( in  the  system  of  coordinates  in  which 
the  shock  wave  moves)  The  frequency  conversion  can  be  quite  large 
due  to  a  change  to  perturbations  of  a  type  different  than  an  incident 
'wave.  Measurements  of  the  frequency  shift  theoretically  allow  us  to 

•  ‘  i 

restore  the  shock  wave  parameters.  ,t 

General  expressions  are  obtained  for  the  coefficients  of  reflec¬ 
tion  from  shock  waves  and  for  the  passage  through  them.  The  final 
formulas  are  given  in  the  simplest  case  of  normal  incidence  to  a 
perpendicular  shock  wave* 

Equating  the  determinant  of  the  system  for  amplitudes  to  zero,- 
we  obtain  a  characteristic  equation  which  should  be  the  basis  of  the 
investigation  of  stability  relative  to  small  wavelik4  perturbations 
of  the  surface  of  discontinuity.  The  region  of  spontaneous  radiation 
of  sound  (resonance  region)  corresponds  to  real  values  ofifrequency 
and  of  wave  vectors.  The  characteristic  equation  for  fast  and  slow 
shock  waves  and  also  for  rotational  discontinuities  have  a  different 
fora. 
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;  FIPSl  l.'  E  OF  TEXT 
1 

i  The  wave  vector 

».  1 

t  r  _ 


2.  Designations 


k  -  -  (veot«,  q), 


The  normal,  to  the  undisturbed,  surface  of  discontinuity  n  *  ("1,0,0) 
The  projection  of  the  wave  veptor  onto  the  plane  of  discontinuity 


q  *  (0,  0  s!h  p,  ?  to*  P). 


The  angle  of  incidence^ 

FIFfT  U'JS  OF  T  ' 


L  L 

art-  to*  -j2-  m  arc  «t_2. . 


The  angle  between  the  direction  of  the  wave  vector  and  the 
magnetic  field  Is  0  .  . 

The  frequency  in  the  frame  of  reference  1,  where  the  surface  of 
the  shock  wave  is  at  rest  is  o>,  whereas  in  the  frame  of  reference 
where  the  fluid  is  at  rest,  it  is  .  The  moduli  of  phase  velocity 

t 

of  perturbation  in  the  direction  x  for  fast  and  slow  magneto-acoustic 
waves  and  Alfven  waves  [4]  ..are  respectively: 


»'+  to- 1 -  to.  I* M- 


(*) 


The  curve  expressing  the  dependence  of  phase  velocity  On  angle 
of  incidence  in  the  plane  of  incidence,  i.e.,  for  P  *  const,  is 


r+to.  i'xto 


(5) 


(in  polar  coordinates) .  I 

The  half- space  in  front  of  the  shock  wave  is  II;  behind  the  shock, 
wave,  it  is  n.  The  magnitudes  pertaining  to  H  are  designated  by  "5, 


and  those  to  II  by  A.  The  Jump  in  magnitude  at  the  discontinuity  is 

(6 ) 


[A)  -.4-4. 


STOP  HESS 
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The  angle  between  the  total  velocity  of  perturbation  +  V  .  H+rxn 
(with  regard,  to  drift)  and  the  velocity  of  this  perturbation  +V  ;  h 
relative  to  the  fluid  at  rest  is  designated  by 


y 

where 


--*■9  •  v  —  r,f6«  «  ±rKM 
•  *r«  r**U« 


t  -fluid  velocity} 
s  « speed  of  sound} 


p  * pressure; 

H  =  magnetic  fluid; 

n  ”•  ■ 

u  * 

p  ‘■density; 
a  *entropy; 
w  =heat  function. 

The  group  velocity  of  perturbation  ( in  the  frame  of  reference 
where  the  discontinuity  surface  is  at  rest) 


r«± 


COS  I  — 


dV 

dm 


(8) 


where  and  are  the  angle  values  at  which  =0  (on  the  Vf  curve) 
fio  and  fi*  are  the  frequencies  radiated  by  the  generator  and 
received  by  the  receiver  which  is  fixed  relative  to  II  (relative  to 
the  fluid  in  front  of  the  shock  wave) ;  the  frequency,  drift  is 

(9) 


The  expression  for  the  perturbed  discontinuity  is 

nexpifor- «i).  (10) 

The  system  of  equations  for  amplitudes  of  divergent  waves  and  of 
the  oscillation  amplitude  of  discontinuity  Tj(iTj  »  6*c/P^)  1® 


-4- 


2 t, 4 1 -  2 u**A r-  7. 


»•! 


(u) 


•  The  index  J.  indicates  the  number  of  the  boundary  equation. 
Coefficients  Tj^  and  Uj^  are  found  by  expanding  the  arbitrary  pertur¬ 
bation  to  the  sum  of  fundamental  waves  :  . 


-  •$,  *  r±M  +  «,« 

•  II  =  b£  i/>±(U  +  4,u  tflfm  • 


(12) 


Here  6p+  is  the  pressure  amplitude  in  the  Vf  wave,  6H  is  the 
z-component  of  the  magnetic  field  in  the  Alfven  wave.  The  index  (l) 
has  the  values  of  1  and  2  and  indicates  whether  the  given  wave  is  with 
the  current  (1)  or  against  the  current  (2)  (if  we  change  to  a  normal 
incidence).  In  (12)  is  the  summation  with  respect  to  +  and  (i)s 


.  _  ■  w»  wjk  —  **(lfu)B 
p  **  **  «-* — i*  “)*  ' 


u* 

“5 


_ i4II~(kII)k 

~ uj— (ku)» 

|  lb.  HI  . 

/4vp  *  Tfc.  HI,  ’  d 


(13) 


lb.  Ill 

ncmr 


If  the  boundary  equations  on  the  discontinuity  are  reduced  to  the 


fora 


Tf> •  i ij  +  (7 j>  +  Pj  +  +Qjlll)3®. 


(14) 


then  the  coefficients  of  equation  ( ll)  will  be  expressed  by  the  combina¬ 
tion  of  magnitudes 


: .  =Pj+  ; 

From  the  continuity  of  mass  flow  ( J  «  l) 


(15) 


ri.  =  (p(«  — q  v)>,  f„-rr«. 
A*=r,/«*.  R,  =  p«.  0,-ft 


(16) 
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From  the  continuity  of  energy  flow  ( J  =  2) 

T f  *=  |p(y  +  »  +  (»  -  q  ')  +  9 *;«  +  ^H-| : 

(4  +  -)  +  pH 


Tat  —  fi 


1  + 


t(t+-) 


(IT) 


R*  =  pr'v7J&  H  +  p(t  +  "’  +  ^),,: 


Q,  =  -t£-  H  —  v - , - 

*  ht 

The  continuity  of  momentum  flow  ( j  *  3>^>5)  yields 

1P>»  =  -2fr,q.{v)1}'s  =  lJ,;  . 

^-1  +  ^M  *-2,r„Vq«. 


(»  H) 


>  -  {p  'i  *•  +  9  r,  («■» - q  •  V)  +  -^-(q  H)J,  Ta  =>r, r, ; 


(19) 


P*=pp-.  n.=  (pr„  Pr,,0).  Q,=  j_A. ]j— v oj  ; 

r«.=  {piA'«+  pr,(«~q-v)  +-^(q  •  H)|,  !■„= 

**•  ~  “V1”’  n»  =  (pp*>  °*  Pr*)'  <**  =  (-‘5r*0>  — 


(18) 

(19) 

(20) 


From  the  continuity  of  tangential  components  of  the  electrical 
field  ( j  -  6,7) 

P*»  =  {  //„  (  u  —  q  .  v)4-  r„  (q  -  H)) ,  r,,  =  0; 


/>.  =  °.  0).  Q.=  (— r, .  r„.  0); 

r;.-(//,(„_,.v)+r,(,H)),  Tn  =  0; 

-  o.  R,  =  <//,,  0,  -//,).  Q,  .  ( — r,,  0.  rMy 


(21) 

(22) 


The  eighth  equation  ( J  =  8,  the  continuity  of  the  normal  compo¬ 
nent  of  the  magnetic  field)  upon  perturbations  of  type  ( 10)  follows 
from  the  two  preceding  equations. 

3-  Results 

We  will  present  a  number  of  results  from  this  examination.^ _ 

-  ■  .  t 

The  law  of  reflection  has  the  form  -  - 


(23) 
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According. to  expression  (7),  the  ends  of  the  phase  velocity 
vectors  for  all  reflected  and  incident  waves  satisfying  expression 
(23),  lie  on  one  circumference  (fr-circumference)  passing  through  the  I 


Pig.  1.  Plotting  of  phase  velocity 
vectors  of  reflected  waves  from  a 
rapid  inclined  shock  wave 
(«/» coiB— for  the 

medium  behind  the  shock  wave.  The 
arrow  indicates  the  points  of  inter¬ 
section  of  the  ^-circumference  with 
curves  V( a)  .  These  points  are  the 
ends  of  the  phase  velocity  vectors  of 
the  incident  waves  (Jl<o  ,  thick 

arrow)  and  of  reflected  wave  (^>  0) 

For  slow  shock  waves  and  magnetohydrody- 
namic  waves  and  for  media  in  front  of 
a  shock  wave  a  similar  plotting  (♦- 
circumference)  is  necessary. 


This  enables  us,  by  means  of  a  cross  section  of  surfaces  V  +  (w) 

r*  ■'  4  i 

and  V  —  (>t)  with  a  plane  of  incidence,  to  construct  the  reflected 

i-  '  •  I 

waves  with  respect  to  a  given  shock  wave;  the  intersection  of  ifr-clr-  : 

cumference  with  V(o)  curves  will  yield  the  ends  of  phase  velocity 


vectors  of  the  unknown  waves.  ,  I 

- 1  2 — - — | 

_ Differentiating  expression  (7)»  we  obtain:  , — — j 

(  ,  ■  1 


whence  it  is  obvious  that  the  separation  into  Incident  and  divergent 
waves  can  be  done  with  respect  to  4he  symbol, Which  is  determined 

from  the  same  diappmw,' ’;.C.  i , Vrfi ' :  \'.y 

The  law  of  refraction  la 


(r«oot}J+A{»}*a  \  (25) 

If  q  •  (▼)  »  Qt .  then  expression  (25)  becomes  "the  law  of  tangents" 


'  tan  }JUn^»  iy)r(.  (26) 

•  .  ! 

In  order  to  find  all  of  the  divergent  waves,  we  must  plot  the  ¥  and 

^-circumferences ,  where  f  and  ¥  are  related  by  condition  ( 25)  ,  and 

must  separate  all  waves  diverging  from  discontinuity  with  respect  to 

the  symbol  .  For  fast  waveB  divergent  waves  do  not  occur  in  H  and 

it  is  sufficient  to  plot  only  the  ^-circumference.* 

•  , 
When  v^-*  0,  ifr-circumference  at  the  limit  occupies  a  position 

tangential  to  the  normal  (x-axis)  at  the  origin  of  coordinates. 

Blotting  of  the  reflected  and  refracted  waves  is  as  before.  Since 

when  vx~*  0  formally  f -*  0,  then  in  place  of  the  derivative  <ty/da 

we  can  introduce  a  derivative  (having  the  opposite  sign)  from  the 

diameter  of  the  periphery  ^  and,  by  this  value,  classify  the  waves  as 

incident  and  divergent. 

The  transition  to  an  Incompressible  fluid  results  In  that 
1%  ~  a  "yr -*0  .  For  this,  however,  there  remains  in  I  the 

♦These  problems  and  also  the  occurrence  of  surface  waves  on 
passage  through  and  reflection  from  a  shock  wave  and  the  role  of 
phase  and  group  velocities  have  been  examined  by  us  in  detail  (1). 
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solution  k(  +  iq,q)  ,  and  in  1  ■  k(—  iq,  q)  ,  representing  a  surface  wave 
Therefore  (ku)  is  finite  when  s->.  *», 


r  1»+  ->  _  Jk  H)k 

»“S  ‘ 


(27) 


The  Expression  U)i  _  t0I«  q  siu«  q  and  the  slow  sound  at  the  limit 

gives  the  Alfven  baric  wave 


77*--* 


...  ■  K'H 

fU  tiiu‘6 


(28) 


The  formulas  for  frequency  drifts  (9)  owing  to  the  Doppler  effect 
and  as  a  result  of  the  conversion  of  some  forms  of  perturbations  to 
other  forms  for  normal  incidence  to  shock  waves  were  presented  in  our 
article  [l],  The  frequency  measurements  can  be  used  for  the  study  of 
shock  waves.  On  the  other  hand,  the  effect  of  change  of  frequency 
with  simultaneous  (generally  speaking  not  small)  change  of  amplitudes 
both  toward  a  decrease  and  toward  an  increase  is  a  new  mechanism  of 
frequency  conversion  with  simultaneous  amplification  ( attenuation)  of 
the  signal  .  This  pertains,  of  course,  to  low  frequencies,  where 
magnetohydrodynamics  is  applicable. 

We  will  introduce  the  coefficients  of  passage  for  a  normal  inci¬ 
dence  of  a  signal  to  a  perpendicular  shock  wave.  For  the  fast  shock 
wave  we  have: 

Ti*  *  *  *  +  Tjt  j.  +  r,.  +  Tj, +  rit  ip,  + 

+  Tji  8 -f  8 // ,2i  =  y  Uui  +'■»; 

=  (7)1 .  Tjt  =  I  >  I, .  r«  =  (j)|t  7j»  —  (j), ; 

Ujk  —  Tih,  i  =  I,  ...  It, 


The  coefficients  of  passage  for  incidence  of  the  V  +  wave 
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(29) 


ip  i  ni  _M_  +  Y  M  —  lx*  (.1/  4. 

ip . v C"  +  A)  '  ‘  - /  v (.'»/TaT; 


„  .»  ^  4-  A1  .V  {.V  +  Y)  -  XIX /  +  A1) 
*/>!.•!  ■  2.M*A'  ‘  S  -/(.I/  -f- .V) A"  h. 


(3P) 


For  incidence  of  the  V  wave 


*P-di  -»/-V(.1/  +  A)  1  -/(.!/+ A). V  .  • 


(31) 


Here  j/  =  f  •/ ,,  a'  =  r ,,  /,,  .1  =  «/*.  l  -M*  ~  A-*, 

/=pT/rr{<-)  =  rJ,/?{/.).  f  =  P  =  (~^. 


On  changing  to  ordinary  hydrodynamics  (N-»  1,  A~*0)  we  obtain 
from  expression  (29)  the  coefficient  of  sound  passage  through  a  shock 
wave  [5]. 
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